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It is well known that for two qubits the upper bounds of the relative entropy of entanglement 
(REE) for a given concurrence as well as the negativity for a given concurrence are reached by pure 
states. We show that, by contrast, there are two-qubit mixed states for which the REE for some 
range of a fixed negativity is higher than that for pure states. Moreover, we demonstrate that a 
mixture of a pure entangled state and pure separable state orthogonal to it is likely to give the 
maximal REE. By noting that the negativity is a measure of entanglement cost under operations 
preserving positivity of partial transpose, our results provide an explicit example of operations such 
that, even though the entanglement cost for an exact preparation is the same, the entanglement 
of distillation of a mixed state can exceed that of pure states. This means that the entanglement 
manipulation via a pure state can result in a larger entanglement loss than that via a mixed state. 
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I. INTRODUCTION 

In quantifying quantum entanglement of two-qubit 
mixed states, various measures are commonly applied 
the relative entropy of entanglement (REE) [3] - 
a measure of the "distance" (or distinguishability) of an 
entangled state from the set of disentangled states, the 
(logarithmic) negativity [1, 3 - a measure of entangle- 
ment cost under operations preserving the positivity of 
partial transpose (PPT) 0], and the concurrence - a 
measure of the entanglement of formation 7] . 

It can be shown analytically that the upper bounds of 
the REE for a given concurrence [1] and of the negativity 
for a given concurrence Q are reached by pure states. 
So one could conjecture that pure states have also the 
highest REE for a given negativity. However, we will 
demonstrate that there are mixed states exhibiting the 
REE for a given negativity (in some range) higher than 
for pure states. Before going into details let us briefly 
describe the entanglement measures. 

II. ENTANGLEMENT MEASURES 

The relative entropy of entanglement in two-qubit sys- 
tems according to Vedral et al. can be defined as 0, Q 

Ejiip) = mm^,evS{p\ \a') - Sip\ \a), (1) 

where the minimum is taken over the set 2? of all sepa- 
rable states (7, and S is the quantum relative entropy 

S{p\\a) = TT{plog^p-p\og^a) (2) 

between states p and a. The REE measures a quasidis- 
tance, say D{p\\a), of the entangled state p from the 



closest separable state (CSS) a. Based on the quantum 
version of Sanov's theorem, one can also interpret the 
REE as a measure of statistical distinguishability of p. 
The choice of 5(^1 |(t) as a candidate for D{p\\a) is by 
no means unique, although this is, to our knowledge, the 
only proposal that coincides for pure states with the von 
Neumann entropy of the reduced density operator. Also 
note that S'(p||(t) is not symmetric and nor does it satisfy 
the triangle inequality; thus it is not a true metric. 

The negativity N(p) for a two-qubit state p can be 
defined by 0, Ell, El : 

Ar(p) =max{0,-2^,„i„}, (3) 

where /imin = mincig(/9'") is the minimal eigenvalue of 
the partial transpose, denoted by F, of p. The negativity 
is directly related to the Peres-Horodecki criterion [3,, ^ . 
The logarithmic negativity, given by \og2[N{p) -1-1], is 
a measure of the entanglement cost Ec{p) under PPT 
operations @, E3l ■ The negativity and logarithmic nega- 
tivity are monotonically related, reaching unity for Bell 
states and vanishing for separable states. So for simplic- 
ity of our further analysis, we use the negativity instead 
of the logarithmic negativity. 

Another measure of entanglement is the entanglement 
of formation Ep{p) [3] or, equivalently for two qubit 
states, the Wootters concurrence Q defined as C{p) = 
max{0, 2maxj \j — A^ }, where the Aj 's stand for the 
square roots of the eigenvalues of p{ay ® ay)p*{py ay), 
and ay is the Pauli spin matrix. 

In the last section, we also analyze the entanglement 
of distillation, Ed{p) Q, a measure of the entanglement 
as the fraction of Bell states that can be distilled using 
the optimal purification protocol. 
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FIG. 1: REE Er{p) as a function of negativity N{p) 
for pure states p^^^ (thick solid curves), the (standard) 
Horodecki states p'^' (dashed curves), the optimal general- 
ized Horodecki states pC-"'^' (uppermost thin solid curves), 
and the Bell-diagonal states p^^^^ (lowest solid curves). 
Marked regions correspond to states exceeding the pure-state 
REE: Blue (red) region shows states ' (p''''^') for which 
E'P>eV (i?r'>max{<',4"'}). 

III. REE WITH FIXED TV FOR PURE AND 
MIXED STATES 

The REE and the entanglement of formation coin- 
cide for pure states, but for mixed states the inequality 
Ef{p) > Eft{p) holds As the concurrence is mono- 
tonically related to the entanglement of formation for an 
arbitrary state, the upper bound of the REE for a given 
concurrence is reached for pure states. On the other 
hand, as shown by Verstraete et al. 0], the negativity 
N{p) of an arbitrary state can never exceed its concur- 
rence C{p). The upper bound of the negativity for a 
given concurrence, i.e., C{p) — N{p), is reached for a 
class of states for which the eigenvector of the partially 
transposed p corresponding to the negative eigenvalue 
is a Bell state d, [ij, 111]. Pure states and also some 
mixed states (including Bell-diagonal states) belong to 
this class. Thus, we see that mixed states cannot give 
higher values of the REE and negativity for a given con- 
currence than those for pure states. In the following we 
will show that the mixed-state REE can exceed the pure- 
state REE for a given negativity. 

An arbitrary two-qubit pure state can be changed by 
local rotations into a state of the form (0 < P < 1) 

l^p) = \/p|oi) + VT~P|io), (4) 

as can be shown by applying the Schmidt decomposi- 



tion [l^. The negativity of \tpp) is simply described by 
N{\%Pp)) = 2y^P(l-P), while the REE, being equal to 
the entanglement of formation, can be given as a function 

of N = NUtPp)) as 

E^PiN)^En{\M)-H,(^l[l + VT^]y (5) 

where iJ2(a^) = ^a;log2 a:^~(l~2;) log2(l — a;) is the binary 
entropy. Equation ([5]) corresponds to the well-known 
Wootters relation between the concurrence and the en- 
tanglement of formation since Ndipp)) = Cdtpp)) 
and Er{\^p)) = EfH^p)). 

In comparison with pure states, let us analyze a mix- 
ture of a maximally entangled state, say the "triplet" 
state 1?/'+) = (|01) + \10))/^/2, and a separable state or- 
thogonal to it, say |00), i.e., [rn |. 

=p|^+)(^+| + (l-p)|00)(00|, (6) 

where the parameter p e (0,1). For brevity, we shall 
refer to ([6]) as the Horodecki state, although alternatively 
it could be named after others (see, e.g., @, i). The 
negativity of the Horodecki state reads as 

iV(pW)- V/(1-P)'+P'-(1-P), (7) 
while the REE as a function of TV = 7V(p(H)) 

can be given 

by Vedral-Plcnio's formula 

(TV) ^ Ep{p^^^) = 2H,{l~p/2)~H,{p)-p 

= ip-2) log2 (1 - p/2) + {l-p) log2(l - m) 

where p = ^2A^(1 -I- N) — N. By comparing the REEs 
for the Horodecki and pure states we observe that 

E'j^\N) > E'j^\N) for < TV < Ny, (9a) 
Sjf ^ (N) < e'^P (TV) for TVy < TV < 1, (9b) 

where Ny = 0.3770... and EP\Ny) = E^ (Ny) = 
0.2279 .. . as shown in the inset plot of Fig. 1. The in- 
equality (j9ap can also be shown by expanding ^ and 
([8]) in power series of N close to zero, then one gets 
eP{N) = TV(1 - y]v72)/ln4 + 0{N^) > and neg- 
ligibly small eP{N) ~ 0{N'^). To show inequality ^ 
more clearly, we can expand ([5]) and ([8]) around — 1 — e 

close to one, then we have E^ (N) = 1 - e/ ln2 + O(e^), 

which is greater than E}^ (N) = 1 — e(l — lne)/ln4 -|- 
©(e^). Thus, a comparison of ([5]) and ([S]) demonstrates 
the main point of this paper: There are mixed states 
having the REE for a given negativity (in some range) 
higher than that of pure states. 

So far, we have analyzed the Horodecki states, which 
correspond to the broken curve in Figs. 1 and 2. Now, 
we will give analytical examples of mixed states more 
entangled than pure states, which could correspond to 
any point in the dark region in Fig. 2(a). The desired 
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FIG. 2: Ranges of values of the REE with a given negativity 
for the following states: (a) p'^ -* given by (|10|l , (b) p'^ ■* given 
by (mi), (c) p'-°^'> given by and (d) p'^^'' given by pi)l . 
Curves are the same as in Fig. 1. 



A^). For p = po = y^2N{l + N) - N, the state 
(fTO| goes into the Horodecki state, given by ([6]). States 
corresponding to all points in the blue region in Fig. 1 
can be generated from p*^^ ^ {p, N) by changing N from 
to Ny and slightly increasing p from the value of pq. By 
choosing properly iV e (0,1) and p e (po,l), the state 
p*-^ ^ {p, N) corresponding to any point in the entire dark 
region in Fig. 2(a) can be generated. It is seen that 
pure and mixed states having the negativity higher than 
that p^^-*, which correspond to the white region under 
the thick solid curve in Fig. 2(a), are not included in 
the family of states p^^ ^ By contrast, dark region in 
Fig. 2(b) corresponds to states that can be obtained 
from pure states {ipp) by mixing them with the separable 
states cr(P) = P|01)(01| + (1 - P)|10)(10| closest to p^^). 
They can be given, in analogy to pO)) . as 



p(P')(P,iV) = (l-:E)|^p)(V;p|+xa(P), 



(14) 



where x = I- 7V/[2^P(1 - P)] for N e (0, 1) and P G 
{P-,P+) with P± = i(l ± Vl - N^)- The bounds on 
P are obtained from the requirement that p^^ ^ (P, N) 
should be a positive semidefinite operator. In special 
cases for P = P±, the mixed state p^^ ^ (P± , N) becomes 
the pure state p(^^(iV). In analogy to the state (fTU)l . 
the Vedral-Plenio theorem guarantees that the CSS for 
p(P )(p^ TV) is the same as for the pure state \ipp). Thus, 
we can calculate the REE for (fT4| arriving at 



mixed states can be generated from the Horodecki state 
p(H) by mixing it with the separable state ct'^' closest to 
p(^^ as follows: 



j^^'Hp,N) = (l-a;)pW +a;(7("\ 



(10) 



defined for N e (0, 1) and p G (y^2N{l + N) - N,l), 
where 



_ (iV + j3)^-2A^(l+iV) 
^~ P^l + N) 

and the corresponding CSS is given by {q — p/2) 

1 

^^(p) = q{l-q) Y,i-iy-''\jA-j){k,l~k\ 
j.k=0 



(11) 



+ (l-g)2|00)(00| + g2|ii)(ii| 



(12) 



By virtue of the Vedral-Plenio theorem [3], state (|12p is 
the CSS for p(H') fo^. ^ny x G (0, 1). Thus, we find that 
the REE for p^^'^P, N) is given by 



£;p(p(H'))^i?)^"^(p,iv) 

q^x log2 X + 2qyi logj ' 



(H')/ 



(13) 



1-9 



y2 iog2 



2/2 



(l-# 



where yi = 1 — qx, t/2 = 1 — 2(7 + q'^x. The choice of 
x, given by JH]), implies that N is just the negativity of 



Ep{P,N) 



(15) 



2P 



where y± = 1 ± x/T~2z and .z = 2P(1 - P)a;(2 - x) = 
2P(1 - P) - iVV2. The REE, given by (US]), for P P± 
goes into ([5]) as expected. The lower bound of the REEs 
for both p(^ ) (P, TV) and ^ (p, N) is the same and given 

by 



E 



(BD) 



(N) ^ Ep{l/2,N) = E),-'{1,N) 

-^-H.C-±^). (16) 



,(H') 



With the help of the Vedral et al. results 0, we can 
conclude that the REE, given by ([T6|) . is the same as for 
a Bell-diagonal state defined by 



p(BD)=^A,|A)(A|, 



(17) 



i=0 



1, maxj- \j 



where \[3i) are the Bell states, 
(1 TV)/2 > 1/2, and N is the negativity N{p^^^^). 
Specifically, the states (HH]) for p = 1 and dH]) for p = 1/2 
go into the following Bell-diagonal states 



p(H')(i,7V) ^ l_Z(|^o)(/3o| + |/32)(/32|) + - 



TV 



l/?3)(/?3|, 
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p(P')(l/2,iV) = ^-±J:L\p,){p,\ + i_i^|/33)(/33|, (18) 

respectively, where the BeU states are given in the fol- 
lowing order |/?2,+fe) = [|0,A:) + - k)]/^. 

The Horodecki state is more entangled than the pure 
state at most at TV = TV' = 0.1539... for which 
maxjv[i^lf^(iV) - e'^P {N)] = 0.0391 . . .. The question 
arises about the highest possible difference for an arbi- 
trary mixed state. This problem is strictly related to 
finding the upper bound of the REE vs negativity. 

IV. NUMERICAL SIMULATIONS 

There has been a long-standing open problem of find- 
ing a closed analytical formula for the REE for two 
qubits, which cor resp onds to finding the a for a given 
entangled state p [l^ and it is argued that the analyti- 
cal solution does not exist [l^. Moreover, there has not 
yet been an efficient numerical method proposed to calcu- 
late the REE for an arbitrary entangled mixed state even 
in case of two qubits. Analytical formulas for the REE 
are known only for some special sets of states with high 
symmetry 0, S, H HO, III, ill, [H . Thus, usually, numer- 
ical methods for calculating the REE have to be applied 
[1, [13, [2^. The complexity of the two-qubit problem 
can be explained by virtue of Caratheodory's theorem, 
which implies that minimalization of the quantum rel- 
ative entropy S{p\\a) should be performed over 79 real 
parameters describing decomposed a [1] . Usually [1, [13] , 
gradient-type algorithms are applied to perform the min- 
imalization. Rehacek and Hradil [23| proposed a method 
resembling a state reconstruction based on the maximum 
likelihood principle. Doherty et al. [25] designed a hier- 
archy of more and more complex operational separability 
criteria for which convex optimization methods (known 
as semidefinite programs) can be applied efficiently. 

Nevertheless, there is a compact-form solution to the 
inverse problem, which characterizes an entangled state 
p for a given full-rank a [H, [2^ : 

p = a-xG{a), (19) 

and 

r -/ for 7. =7, . . 

= < 7.-7, for-T-T^-Y- ^ ' 

and \(f>) is the kernel of , while \i > and ji are eigen- 
states and eigenvalues of a, respectively. Thus, the REE 
is given by 

Enip) = S{<j) - S{p) + xtT[{\^){cj,\ fa\og, a] , (22) 

where S{-) is the von Neumann entropy. In the following 
a^max denotes the largest x in (|19p , for which p is a positive 
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FIG. 3: REE Er{p) versus negativity N{p) for randomly- 
generated states p of rank 2 (upper figure) and rank 3 (lower 
figure). The upper- and lowermost curves correspond to the 
optimal generalized Horodecki states and Bell-diagonal states, 
respectively. Random rank-4 states lie in the same range. 



operator. The solution can be applied also for lower- 
rank CSSs in a limiting sequence from a full-rank state 
by noting that the REE is a continuous function. 

We have applied two approaches in our numerical sim- 
ulations. In the standard approach, the states are chosen 
at random and their and N are calculated numerically 
using a simplex search method without using numerical 
or analytic gradients. However, given the fact that no 
closed formula exists for Er [l^, [l3| , it is more logical to 
choose random states (call them a) on the boundary of 
the separability region and generate entangled states p 
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for which a is the CSS by applying Eqs. ([T9l)-(l22l). The 
latter method is faster by three orders of magnitude than 
the standard one. Figure 3 shows the results of our simu- 
lations for random rank-2 and rank-3 states. For brevity, 
we have omitted a similar figure for random rank-4 states. 
The simulations confirm our analytical predictions that 
the mixed-state REE can exceed the pure-state REE but 
also indicate lower and upper bounds of the REE vs neg- 
ativity. 



V. REE FOR THE GENERALIZED 
HORODECKI STATES 

Our numerical simulations, partially shown in Fig. 3, 
indicate that the upper bound Eji{N) can be given by 
the rank-2 generalized Horodecki state defined as 

follows fii: 



JGH) _ 



p|V'p)(^p| + (l-p)|00)(00|, (23) 



where \'4'p) is given by ffl and p,P £ (Oil)- In the 
special case of P — 1/2, p(C!H) j-g^jucgg iq the standard 
Horodecki state, while for p = 1 corresponds to a 

pure state \ipp). Note that the negativity A'' = iV(p(°^') 
is simply given by: 



N=^{1- pY + VP(l _ F) - (1 - p). (24) 
By the inversion of this equation, 

' (25) 



P = 



1 

Yp 



± _ _ 2N{l-p) 



one can express the state p('^H) = p^*^^) (p, iV), given by 
(|23p . as a function of negativity N and parameter p > 
Po{N) = y^2N{l + N)-N. The state p'-^^^ is a special 
case of a more general state 



'n ON 

r2 y 

y r3 

.0 Oy 



(26) 



for which the CSS is given by: 




(27) 



where Y — y/RiRl. Clearly, by assuming y — ^JT^r^^ 
the state p'^-' is reduced into the generalized Horodecki 
state p(^H) with n = 1 — p, r2 = Pp and ra = (1 — P)p- 
States p'-^-' and a'^-' are related by the following relations 
assuming for convenience that R\ > R4 [l8l |: 



r2 = R2 



i?2i?3 + 2F')- 



2i?4 



(i?2-i?3) (28) 



together with ri = i?i — i?4, = 1 — ri — r2, and y = 
-[{r2-R2){R2-R3)-2iRi+R2) Ri]/i2Y) given in te rms 
of the auxiliary functions z ~ y/ {R2 — Rs)^ + 4Y^ and 
L = \ii{R2 + R3-z)-ln{R2+R3 + z). Moreover, Xmax = 
(i?i + i?4)/i?i if the condition y = y^rJFa is satisfied for 
a given choice of {i?^}. 

These equations can easily be inverted for P — 1/2, 
which leads to the solution given by (|12p for the standard 
Horodecki state p^^\ By contrast, due to presence of 
logarithmic functions of nonlinear combinations {Ri} in 
the equations for {ri}, it looks impossible to invert the 
equations in order to express all {Ri} in terms of {ri} for 
the generalized Horodecki state if P ^ 0, ^, 1. Thus, we 
can only give a formula for the REE for p*-*^^^ with {r^} 
as a function of {Ri}: 

= -H2{ri) - n log2 Ri - fl log2 A_ - fl log2 A+, 

(29) 

where /± = J\f±[{X± - R3)V^+Y^, A± = i(i?2 + 
i?3 ± z) and AA± [(A± - P3)' + r']-^/'. 

In any case, a multivariable numerical procedure for 
finding the CSS cr*^^^-' can be reduced to a single- variable 
problem; namely, we can express Ri (for i = 2,3,4) in 
terms of r\ and R\ as follows: 



1 



Ri 




R4 


= i?i - ri. 


R3 






i^3 



2r2 - 4Ri - VS), 



(30) 



where 



r2,3 



1 - ri T Vil-rif - N{N + 2ri) 



S = (3ri + 1)2 -4r2r3 -8i?i(ri -f' 1) 



-16V^i(i?i-ri) 



r2r3- 



(31) 



Thus, to completely determine (t'-'^^' for a given p^'^^\ 
it is enough to find i?i (> ri) by numerically solving the 
single- variable equation ([28]) with all the other variables 
defined above. 

A related problem is to find the optimal generalized 
Horodecki state pC-"^^), defined as p^*^^) for a given N 
and such p, denoted by Popt{N), for which the REE is 
maximized: 

4°™) (TV) ^ £;;^[p(o™)(7V)] 

^ En[p^''''HPopt{N),N)] 

= max £;;^[p(™)(p(7V),7V)]. (32) 

P>Po{N) 

The parameter popt can be found numerically by the pro- 
cedure described above. On the other hand, we have 
found a fairly good approximation of Popt for < A'^ < 
0.527 given by 



PoptiN) 



-N - 



7 

11' 



(33) 
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such that ER[p'^'^^\popt{N),N)] deviates by the order 

10~^ from the precise value of e'^°'^^\n). We also 
find that popt becomes 1, so the optimum generalized 
Horodecki state becomes a pure state |V'p) for iV>0.53. 
It is worth noting that the precision of our numeri- 
cal calculations of the REE is 10"^" 10"**, and 
ma.x{E'^^\E'^P) is smaller than up to 0.0148 (at 

N = 0.377), so it can clearly be distinguished from the 
numerical noise. 

The REE for the generalized Horodecki states p^*^^) as 
a function of N for arbitrary values of p correspond to 
the dark region in Fig. 2(c). In analogy with the states 
p(H ) ^ given by (fTO|) , one can also define a class of more 
general states by mixing pC^^' with its CSS cr'^'^^^ , given 
by (fT2|l . as follows: 

p(GH')^(l_^yGH)_^^^(GH)^ (34) 

where x G (0, 1). As is seen in Fig. 2(d) in comparison to 
Figs. 3(a) and 3(b), the REE vs for p^^iH') covers to 
whole region of the values for randomly generated states. 

We conjecture that for any two-qubit state p described 
by the REE e\^\n) = Eji(p) as a function of the nega- 
tivity N = N{p), the following inequalities are satisfied: 

^(OGH) ^ ^(p) ^ ^(BD) ^ ^gg^ 

which simplify to 

E^P (N) > E^p [N) > E^^^^ (N) (36) 

for iV>0.53, where E^iN) and E'j^^\N) are given by 

([5]) and (fT6|) . while E^^^\n) is found numerically by 
the described method using Eqs. (P^ - ([?T|) . As a partial 
analytical support of our conjectures, the extremal condi- 
tions for the REE with a fixed N for g^^j p(BD) are 
examined in the next section. We have also performed 
a numerical analysis, as discussed in Sec. VI, to provide 
another support of validity of the conjecture. We have 
generated altogether a few million random states p of a 
fixed rank (2,3, and 4) and calculated the negativity and 
REE for each of them. 



VI. SOME EXTREMAL CONDITIONS FOR 
REE WITH FIXED TV 

In the following, we show analytically that the Bell- 
diagonal states and the generalized Horodecki states, 
thus also pure states and the standard Horodecki states, 
satisfy some extremal conditions for the REE with a fixed 
N implied by a Lagrange-multiplier method. Since neg- 
ativity for a given state p is given by 

Nip) = -2min[(^'|pn^'>] (37) 

= -2min[trp(|V')(V''l)n ^-2Ml^)(V'|)n> 



where \ip) is the optimal state, let us consider the follow- 
ing Lagrange function: 

£ = trplog2p-trplog2a + /(^trp(|V')(V'l)'^ + y) , (38) 

where / is a Lagrange multiplier. For a small deviation 
of 

p -> p + A ~ (trA)p, (39) 

where A is an arbitrary (but small) operator on the sup- 
port space of p [denoted by supp(p) hereafter], we have 

C ^ /:+trA[log2 p-log2 a+m){^\f~ERip)+^-N{p)] . 

(40) 

Since A is an arbitrary operator on supp(p), the following 
extremal condition is obtained: 

P[log2 p - log2 a + imiijlf - Enip) + ^A^(p)]P = 0, 

(41) 

where P is the projector to supp(p). Moreover, the ex- 
tremal condition for £ with respect to {ip) leads to the 
extremal condition for negativity, and thus |^) must be 
the eigenstate corresponding to a negative eigenvalue of 
p^. Therefore, it is found that p, its closest separable 
state tJ, and the eigenstate corresponding to a nega- 
tive eigenvalue of p^ should satisfy Eq. (|^T|l . 

Now let us consider the case where p is a mixed 
state of rank 2, i.e., p — Ai|ei)(ei| -I- A2|e2)(e2|, where 
{Xi} are nonzero eigenvalues of p, and jej) are corre- 
sponding eigenstates. The projector P is then P = 
|ei)(ei| + |e2)(e2|, and as a result the extremal condition 
of Eq. (|^ becomes 

(ei|log2a|e2)=Z(ei|(|V)(^|)^|e2), (42) 

and 

log2 Ai - (eil log2 (7|ei) + Z(ei|(|V')(^|)^|ei) 

= Enip) - i^Nip), (43) 

l0g2 A2 - (62 I l0g2 a 1 62 ) + Z (62 I ( I V') 1)^ 1 62) 

= Enip) - ^Nip). (44) 

However, Eqs. and (I44p are not independent of each 
other. Indeed, for Ai 7^ and A2 7^ 0, these equations 
are equivalent to 

Ailog2Ai - Ai(ei|log2(T|ei) +Ai;(6i|(|V')(V'l)^|ei) 

= XiEnip)-Xi^-Nip), (45) 

A2l0g2 A2 - A2(62| I0g2<7|e2) -|- A2^ (62 | ( | V') ("0 1 )^ 1 62) 

= X2Enip)~X2^-Nip),iA6) 

and it is found that the sum of these equations is au- 
tomatically satisfied. Therefore, the extremal conditions 
for rank-2 states are Eqs. dH]) and (g^l) [or Eqs. (gH) 
and dH])]. 
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A. Bell-diagonal states 

For the rank-2 Bell-diagonal states, 

[^(BD)^^(BD)]^[^(BD)^(|^^^^|)r]^0^ (47) 

and hence (ei| log2 cr|e2) = (ei KIV') ('0l)'"|e2) = holds 
again. Equation is then satisfied for any I, and the 
extremal conditions are satisfied. 

B. Horodecki states 

For the standard Horodecki state, defined by ([6]), with 
its CSS given by (fT2|) . we have: 



(pW)r = g(|V,+)(V;+| + |V^-)(V'-| + |00)(ll| 

+ |ll)(00|) + (l-p)|00)(00|, (48) 



— |oo)-yiT|ii) 



and 



mi^Plf = .s_|00)(00| + ,s+|ll)(ll| 

+ |(l^->(V'-|-|^+)(^+l), (49) 

where 1^ ;=^) = (|01) ± |10))/V2, s± = [l±t{l-p)]/2 and 
t = l/^2p2_2p+l. 

Since (ei|log2a|e2) = (ei|(|^)(V'|)^|e2) = 0, Eq. (gH) 
is satisfied for any I. Since there is only one relation of 
Eq. P5)) for I, the extremal conditions are necessarily 
satisfied. 



C. Generalized Horodecki states 

Here, we show that the generalized Horodecki states 
are extremal. The point is that only two extremal condi- 
tions should be satisfied for the states: (i) given by ([42]) 
and (ii) given by either of Eqs. or, e.g., the 

difference of Eqs. and Condition (ii) is a lin- 

ear function of I, so it can easily be solved for /. The 
question is whether the found / also satisfies condition 
(i) or whether left- and right-hand sides (LHS and RHS) 
of (i) are equal to zero. In the following we show that 
the latter case is satisfied for the generalized Horodecki 
states p^'^^\ and thus also for the optimal states p^^*^^), 
the standard Horodecki states p^^-* , and pure states p^^-* . 

For simplicity, we use the notation of Eqs. and 
(I27p with the condition = r2r^, which guarantees that 
p(^) and CT*^^' become p^*^^) and cr^^^' , respectively. One 
finds that 



IV) = A/-(-ff|00)+2y|ll)), 
m{Mf - AA2[5'|00)(00| + 4y2|li)(li| 

-25y(|01)(10| + |10)(01|)], 



where g = \/r\+ 4y^ — ri and J\f = 1/ y^g^ + 4y^. On 
the other hand, logj a^'^^^ can be calculated through the 
eigenvalue decomposition 



^(GH) 



i?i|00>(00| + i?4|ll)(ll| 
+A+|A+)(A+| + A_|A_)(A_|, (51) 



where 



i?2 + i?3 ± V(^2-i?3)2+4y2 



|A±) = A4[(A±-i?3)|01) + y|10)], 



(52) 



with 7V± = [(A± - i?3)2 + r2]-i/2^ Moreover, for the 
nonzero eigenvalues, the eigenvectors of p^^^) are found 
tobe|ei) = |00) and|e2) = (l/VFT7J)(y|01)+r3|10)). 
Thus, it is seen that both LHS and RHS of condition ((i^ 
are equal to zero. The second condition is satisfied by 
choosing 



log2 ri -hlog2 i?i) 



if + m-n) 



(53) 



(50) 



where / = y^r^ + Ay^. In any case, even without know- 

ing exphcitly and {Ri} in terms of {r^}, we have 

showed that the generalized Horodecki states satisfy the 
extremal conditions. 



VII. CONCLUSIONS 

In conclusion, we have demonstrated that there are 
mixed states that have the REE in some range of a fixed 
negativity higher than the pure-state REE for the same 
negativity. This is somewhat surprising, since mixed 
states can neither exhibit the REE for a given concur- 
rence nor negativity for a given concurrence higher than 
those for pure states. By applying Lagrange multipli- 
ers, we have also shown that the Bell-diagonal states, 
pure states, but also the so-called generalized Horodecki 
states, which are mixtures of a pure entangled state and 
pure separable state orthogonal to it, satisfy some ex- 
tremal conditions for the REE with a fixed negativity. 

Our findings implicitly show another fact. For a given 
negativity, the entanglement of distillation (ED) Ejj of 
the Horodecki state p^^^ can be larger than that of a pure 
state p^^\ For example, N = 0.1 and p = 0.37, one gets 
£:c(p(P)) = EbXp^^^) = 0.025 and Ed{p^^'^) > 0.034. 
Here, the lower bound of the ED is given by p^/4 (via 
the direct method shown in Q). This lower bound can 
be slightly improved by a method discussed in Ref. [27j . 
The point is that the logarithmic negativity is equal to 
a PPT entanglement cost for an exact preparation, the 
REE is equal to a PPT distillable entanglement for pure 
states, and the ED is a lower bound of a PPT distill- 
able entanglement. So our findings provide an explicit 
example of PPT operations where, even though the en- 
tanglement cost for an exact preparation is the same, the 
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ED of a mixed state can exceed that of pure states. In 
other words, the entanglement manipulation via a pure 
state can result in a larger entanglement loss than that 
via a mixed state. 
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